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ABSTRACT 


Studies  of  the  effects  of  sinusoidal  surface  waves  on 
the  propagation  of  sound  in  a  waveguide  reveal  a  significant 
resonance  effect  under  certain  "matching"  conditions  of  the 
wavelength  of  the  surface  wave  and  the  wavelength  of  the 
acoustical  sound  field. 

An  acoustic  waveguide  with  constant  cross-section  and 
pressure  release  boundaries  was  designed  with  a  termination 
to  minimize  reflection  of  both  surface  and  acoustic  waves 
so  that  traveling  surface  waves  and  traveling  acoustic  waves 
could  be  generated  simultaneously. 

With  acoustic  and  surface  waves  present,  the  existence 
of  predicted  resonances  was  verified  and  a  strong  resonance 
was  examined  for  comparison  with  theory.   The  predicted 
effects  have  been  qualitatively  and  quantitatively  supported. 
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I.   INTRODUCTION 

The  propagation  of  acoustic  waves  in  ducts  of  non-uniform 
cross-section  has  been  of  interest  to  theoreticians  for  the 
past  hundred-odd  years.   In  particular,  Lord  Rayleigh  [1] 
investigated  plane-wave  reflection  from  a  sinusoidal  surface 
as  early  as  1875.  J-    C.  Samuels  [2]  presented  theory  des- 
cribing the  propagation  of  harmonic  signals  in  two-dimensional 
ducts  with  slightly  rough  rigid  walls.   R.  F.  Salant  [3,4] 
analytically  investigated  acoustic  plane-wave  propagation 
parallel  to  two  rigid  sinusoidal  walls.   Relationships  be- 
tween surface  wavenumber  and  acoustic  wavenumber  were  des- 
cribed, and  decay  of  the  disturbance  as  a  function  of  distance 
from  the  surface  was  discussed.   A.  H.  Nayfeh  [5]  developed 
a  perturbation  solution  for  acoustic  wave  propagation  in  a 
hard-walled  two-dimensional  duct  whose  walls  had  weak  sinu- 
soidal undulations.   Resonance  conditions  for  certain  rela- 
tionships between  wavenumbers  of  the  wall  undulation  and  the 
acoustic  mode  were  discussed.   W.  E.  Jordan,  Jr.  [6]  compared 
the  frequency  spectra  of  acoustic  fluctuation  and  surface 
fluctuation  created  by  wind-generated  waves. 

More  recent  work  by  R.  H.  Ebert  [7]  examined  the  influence 
of  standing  gravity  waves  on  standing  acoustic  waves  in  a 
water-filled  waveguide  with  pressure  release  boundaries. 
It  was  confirmed  that  an  infinite  number  of  resonances  exist 
when  the  surface  wavelength  is  one-half  the  acoustic  wavelength 
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in  the  x-direction.   Differences  between  experiment  and 
theory  were  attributed  primarily  to  excessive  surface-wave 
amplitude,  non-uniform  standing-wave  patterns,  and  the 
assumption  that  acoustic  and  surface  waves  were  undamped. 
The  present  work  simplifies  the  experimental  problem 
by  examining  the  case  of  a  traveling  acoustic  wave  perturbed 
by  a  traveling  surface  wave  in  a  waveguide  of  uniform  cross- 
section  with  pressure-release  boundaries.   The  investigation 
includes  verification  of  theoretical  resonances  in  higher 
modes  and  study  of  a  strong  resonance  for  qualitative 
analysis. 
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II.   THEORY 

A.   THE  NON-ABSORPTIVE  CASE  -  RESONANCE  CONDITION  [1-8] 

A  rectangular  waveguide  with  cross-sectional  dimensions 

I     and  i     is  excited  by  a  sound  source  at  x  =  0.   The  wave- 

z      y 

guide  is  acoustically  terminated  at  large  x  so  that  there  is 
no  reflection  of  the  incident  acoustic  energy  back  toward 
x  =  0.   All  other  boundaries  are  pressure  release. 

A  well-known  solution  of  the  wave  equation  for  the  ve- 
locity potential  with  these  boundaries  is 


*  =  sin(k  z)sin(k  y)cos(wt  -  k  x) .  (2.1) 

z       y  x 


The   k's   must    satisfy 


k      =  m\/l     and  k     ■  pir/JSi  (2.2) 

z  z  y  y 


where  n,p   =    1,2,3,    ... 

Let  the  upper  boundary  be  perturbed  by  a  traveling  sur- 
face wave  of  amplitude  A.   If  a  perturbation  constant  is 

defined  as  e  =  A/fL     the  boundary  conditions  are 

z 


$  =  0   at  J 


z  =  0 

(2.3) 

Z  -  £   +  eF(x,t) 
z 


where 


F(x,t)  =  -  I   cos(fit  -  yx).  (2.3) 

z 
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FIGURE  1.   THE  PERTURBED  WAVEGUIDE 

Assume  there  exists  a  solution 


$  =  $  +   E   e  $ 

°   n=l    n 


(2.4) 


where  *  .  given  by  Eq.  2.1,  is  the  solution  with  an  unper- 
turbed boundary,  z  =  £_.   The  boundary  conditions  are 

z 


*  =  0  at   - 


z  =  0 
z  =  I. 


Expansion  of  $(z)  about  z  =  £   in  a  Taylor's  expansion  and 
application  of  the  boundary  conditions  leads  to 


94> 
o 

az 


F(x,t) 


(2.5) 


8*. 

Tz" 


F(x,t) 


(2.6) 


1*4 


8$, 


dz 


i    92*u 

F(x,t)  -  ±r  —J-         F^x.t) 

I  3z 

z 


1  9\ 


^-^ 


F3(x,t) 


(2.7) 


It  can  be  seen  that  all  $  have  the  same  y-dependence  as  $ 
The  y-dependent  component  of  <J>  can  be  written  as  Y(k  y) 
where  Y(k  y)  is  the  solution  to  the  Helmholtz  equation 


2-~   Y  +  k   Y  =  0. 

^7       y 


(2.8) 


For  the  non-absorptive  case,  $  and  its  components  must 
satisfy  the  lossless  wave  equation: 


CT$  =  0, 


(2.9) 


n  *n  =  o. 


(2.10) 


Examination   of   the    first-order  approximation   gives 


=    (k   *    )Y rt5-{cos[(u+Ji)t    -    (k  +y)x] 


(2.11) 


+    cos[(u-n)t    -    (k      -   y)x]}. 
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For  this  boundary  condition  an  acceptable  solution  is 

,        sinK  z 

♦.  =  ±(k  a  )Y{ 1 —  cos[(w+n)t  -  (k  +y)x] 

1    d      z  z    sinK  £  x 

z  z 


sinK  z 

+ —  cos[(u-n)t  -  (k  -y)x]} 

sinK~S. 
z  z 


where  the  propagation  constant  in  the  z-direction  is 


(2.12) 


1 
l±2(k  /y)  2 


K-  =  k  [1 ft-,-]   .  (2.13) 

Z     Z       (kz/Y) 

+ 
K   can  be  real  or  imaginary  depending  on  the  relative  values 

E 

2 
of  k  /y   and  (k  /y)  •   There  is  a  condition  of  resonance  when 

4. 

K~  is  real  and 
z 


K~Z      =  m7T        where  m  =  1,2,3  ....        (2.1*1) 

£4  C-i 


Note  that  K   is  imaginary  in  most  of  the  region  where  k  /y   >    0 

and  therefore  there  is  no  resonance  due  to  K  .   K  will  be 

z    z 

referred  to  hereafter  simply  as  K  . 

Li 


Equation  2.2  stated 


k   =  y—  where  n  =  1,2,3  •••• 

z 


therefore  a  resonance  will  occur  where 


2  .      1  -  2(k  /y) 
(kz/Y)2 
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p 

The  relationship  between  k  A  and  (k  /y)      for  the  resonance 

X  Ct 

condition  can  now  be  written 


T  =  "  l[(kz/Y)2  "  1]  +  l(Ff)2(T)2  •       (2-l6) 


Since  resonances  occur  only  when  K   is  real, 


1  -  (2k  /y) 

*5 —  <  1  .  (2.17) 

(kz/Y)^   " 

As  can  be  seen  in  Figure  2  the  theoretical  plot  of  k  /y 

2 
vs  (k _/y)   is  a  linear  relationship  for  the  various 

perturbation  modes. 

B.   THE  ABSORPTIVE  CASE  -  PROPERTIES  NEAR  RESONANCE  [8] 

The  wave  equations  for  the  acoustic  potentials  in  the 
absorptive  case  can  be  written  as 


(°2  -2^)$o  =  0  <2'18> 

(D2  -  2^1  ^r)  •  -  0  (2.19) 

cQ  8t   n 


where  a  and  a  are  the  coefficients  of  absorption  due  to 
on 

all  acoustical  losses,  and  can  be  considered  approximately 
equal  in  magnitude  if  *  and  <!>  have  similar  form. 
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(k2/y)' 

FIGURE    2.       THEORETICAL    RESONANCE    CURVES 
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Let    $  and   <J>     be   of  the   form 
n 


Y(kyy)Z(kzz)exp[i(u)t   -  k_xx)] 


where  k   ■  k  -  in  .   It  can  be  shown  that 
—x    x     'n 


a 
_  0)   n 


a 


n 


cos  a 


I  ■       i        *       -       ;    i 

FIGURE  3.   RELATIONSHIP  OF  k  AND  k 


(2.20) 


(2.21) 


where  n  is  the  acoustic  absorption  coefficient  in  the  x- 
direction  and  can  be  determined  empirically. 

The  classical  solution,  Eq.  2.1,  now  becomes 


*  =  exp(-n  x)Y  (k  y)sink  zexp[iUt  -  k  x)]  .     (2.22) 


If  there  is  surface  wave  attenuation,  Eq.  2.3  is  rewritten 


F(x,t)  =  -I   exp(-6x)cos(Qt  -  yx) 


(2.23) 
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where  the  surface  wave  attenuation  constant  6  can  also  be 
determined  empirically. 

To  determine  the  damping  in  the  z-direction,  assume  $ 
has  the  form 


$  =  Y(k  y)exp(-iK  z)exp[i(a)-ft)t]exp[-i(k  -y)x] 
y  s  x 

(2.211) 
exp[-(n+<S)x] 


where  K   =  K   -IB. 
— z     z 

Substitution  of  this  equation  into  Eq.  2.19  shows  that  the 
spatial  decay  in  the  z-direction,  g,  is  given  by 


ouk  -  (n+6)(k  -  y) 

3  =  — g •  (2.25) 

z 


The  perturbed  solution,  Eq.  2.12  then  becomes 

,       sinK  z 
$1  =  ?(kzJlz)YsinKZ£  Gxp[-(n+6)x]exp{i[(oj-^)t-(kx-y)x]} 

— z  z 

(2.26) 

where 


AK   =  K   -  k    and   A? =  tan"1  Ay  „   .  (2.27) 

z    z    z  AK  a 

z  z 


The  first  order  solution  for  the  acoustic  potential, 

•  ■  #  +  *.  ,  in  the  vicinity  of  a  resonance  near  an  antinode 

in  the  transverse  directions  (y  and  z)  can  now  be  written 

as 

20 


*  = 


where 


Yexp(-nx)exp[i(wt-k  x)]{sinkz+^Rexp(iAf)exp[i(nt-Yx)]} 

(2.28) 


R  = 


siniTz/£ 


(AK  I    )2   +  (BA  )2 
z  z        z 


z  /   0,1, 


(2.29) 


and 


e  =  e  exp(-6x) 


(2.30) 


The  total  solution  can  be  viewed  as  a  sum  of  a  fixed  and 
a  rotating  phasor  as  shown  in  Figure  4.   The  modulation 
parameter,  M,  is  defined  to  be 


M 


=  T5-TTE  Rexp(-6x)  . 

c.        O 


(2.31) 


FIGURE  4.   PHASOR  REPRESENTATION  OF  SOLUTION 
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C.   COMPUTATION  OF  ACOUSTIC  FREQUENCY  AND  SURFACE  WAVE 
NUMBER  -  RESONANCE  CONDITION 

Substitution  of  Eq.  2.2  into  Eq.  2.16  yields 


k/Y  ■  |[m2  -  n2]^)2  +  |  (2.32) 


which  states  the  relationship  between  the  acoustic  wave 
number  in  the  x-direction  and  the  surface  wave  number  for  a 
given  perturbation  ratio,  m/n.   The  relationship  between  k 

A 

and  the  acoustic  angular  frequency  u  is  given  by 


kx  ■  C<t>2  -  \2    -   x/l2  •  <2-33> 


When  %      =    I    j  application  of  Eq.  2.2  gives 

1 

kx  =  [(|)2  -  (n2  +  P2)(JL)2]2  (2.3*0 

z 


which  can  be  written  as 

1 

S  =  Ck  2  +  (n2  +  p2)(^)2]   .  (2.35) 

ex  xz 


The  angular  cutoff  frequency  of  this  square  waveguide  is 


1 

^  =  [(n2  +  P2)(^)2]2  .  (2.36) 

c  *z 
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Squaring  Eq.  2.36  and  substitution  into  Eq.  2.36  yields 

1 

I"  [kx2  +  Hr)2]2  '  «-3T) 


which  can  be  written  as 


X 

w  =  [kx2c2  +  ^c2]2.  (2.38) 


In  terms  of  the  acoustic  frequency,  Eq.  2.38  becomes 

1 
f  =  [(2?)  c2  +  C^f)2]  (2.39) 

1 

2 

=  [(r-)2  +  f.2J  •  ^2-z,°) 

x      c 


Thus  for  a  given  y,  the  two  values  of  k  given  by 

A 

Eq.  2.32  could  result  in  the  same  acoustic  frequency  for 
different  combinations  of  normal  mode,  (n,p)  and  perturbation 
ratio,  m/n.   The  acoustic  frequency  and  surface  wave  number 
may  be  degenerate  at  one  point  only  or  they  may  be  coincident 
throughout  the  range  of  resonant  acoustic  frequency  and 
surface  wave  number.   In  this  latter  case  the  different  com- 
binations of  the  normal  mode  and  perturbation  ratio  are 
termed  degenerate  modes.   For  example,  a  comparison  of  Figure 
1^,  for  the  (1,1)  normal  mode  with  perturbation  ratio 
m/n  =  2/1  with  Figure  16  for  the  (2,1)  normal  mode  with 
perturbation  ratio  m/n  =  p-,  reveals  degenerate  modes. 
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III.   EXPERIMENT 

A.   DESCRIPTION  OF  THE  SYSTEM 

An  acoustic  waveguide  with  pressure-release  boundaries 
was  designed  for  the  experiment  (see  Figures  5,  6,  and  7). 
The  size  was  chosen  to  provide  a  uniform  water  column  10 
centimeters  square  and  1.52  meters  in  length.   These  dimen- 
sions produced  a  calculated  cut-off  frequency  for  the  lowest 
mode  of  10.466  kHz,  for  a  speed  of  sound  in  water  of 
c  =  1^80  m/sec.   The  walls,  constructed  of  one-inch  thick 
polyurethane  foam,  thinly  coated  with  liquid  neoprene  to 
provide  strength  and  minimize  wetting  of  the  foam,  closely 
approximated  pressure  release  surfaces  as  verified  by 
pressure  probe  measurements. 

A  waveguide  termination,  designed  to  minimize  reflections 
of  the  acoustic  wave,  consisted  of  a  20  inch  section  of  an 
exponential  horn  and  a  16  inch  V-shaped  closure.   The  horn 
was  constructed  of  styrofoam  and  the  closure  was  lined  with 
sound-absorbent  aluminum-impregnated  rubber.   An  angle  of 
60°  in  the  closure  assured  that  a  plane  wave  propagating  in 
a  straight  path  down  the  waveguide  would  suffer  at  least 
two  reflections  before  being  reflected  back  into  the  channel. 

To  reduce  surface  wave  reflections,  a  "beach"  was  con- 
structed, consisting  of  a  wire  screen  tray  cut  to  fit  the 
60°  closure,  and  penetrating  the  surface  of  the  water  to  a 
depth  of  about  one  centimeter.   The  tray  was  filled  with 
stainless  steel  ribbon  shavings  to  dissipate  the  surface 
wave  energy. 

2k 


FIGURE  5.   PLAN  VIEW  OF  ACOUSTIC  WAVEGUIDE 
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Constant  Voltage  xfmr 
Reg  ac  Pwr  Supply 
Amp  #1 

Pwr  Amp  #1  and  #2 
Electronic  Filter 

VTVM  #1  and  #2 

CRO  #1 

CRO  #2 

Audio  Osc 

WAVETEK 

Laser 

300/600  Vdc  Bias  Battery 

Variac 

Probe 

Xdcr 

Wave  Maker 


SOLA  Electric  Company  #3097*1 

CENCO  #  79642 

Hewlett  Packard  Model  465A 

Hewlett  Packard  Model  467A 

Spencer-Kennedy  Laboratories,  Inc., 
Model  302 

Hewlett  Packard  Model  400D 

Hewlett  Packard  Model  120B 

Hughes  MEMOSCOPE  Model  105 

CENCO  #80593 

WAVETEK  Multi-purpose  VCG  Model  116 

C.  W.  Radiation,  Inc.,  0.5  mW  HeNe 
Laser 

Burgess  U200  300V  "B"  Battery 

Superior  Electric  Co.  POWERSTAT 

CELESCO  Industries  Hydrophone 
Model  LC5-2 

Locally  Manufactured  MYLAR 
Electrostatic  Source 

Locally  Manufactured  Surface 
Agitator 


FIGURE  7.   LIST  OF  EQUIPMENT  SHOWN  IN  FIGURE  6 
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The  surface  waves  were  generated  by  means  of  a  variable- 
speed  ac  motor  rotating  an  adjustable  cam  attached  to  a 
driving  arm.   The  driving  arm  was  connected  to  a  wedge  of 
plexiglas  which  penetrated  the  surface  of  the  water  and 
generated  surface  waves  by  vertical  oscillatory  motion. 
The  surface  wave  frequency  could  be  adjusted  in  the  range 
2-6  Hz  by  varying  the  supply  voltage  to  the  motor.   The 
frequency  was  measured  by  means  of  a  cam-driven  microswitch 
which  actuated  an  electromechanical  counter.   Wave  amplitude 
could  be  adjusted  in  the  range  0-5  mm  by  varying  the  eccen- 
tricity of  the  cam  and  was  determined  by  a  slope-detecting 
optical  system  described  in  Figure  8.   A  laser  was  used  to 
provide  an  intense,  collimated  beam  of  light  which  was  re- 
flected onto  the  perturbed  surface  of  the  water  by  a  set  of 
mirrors  which  could  be  traversed  in  the  x-direction.   The 
reflected  light  was  projected  to  a  screen  over  a  long  path 
length.   The  system  was  calibrated  by  mounting  a  small  mirror 
on  a  one  meter  bar  and  then  reflecting  the  light  beam  from 
this  mirror  instead  of  the  water  surface.   The  bar  was  then 
set  at  known  angles  while  the  deflection  of  the  beam  on  the 
screen  was  measured.   The  observed  deflection  was  proportional 
to  the  angle  of  inclination  of  the  bar  over  the  range  of 
interest. 

The  acoustic  field  was  generated  by  a  plane  MYLAR 
electrostatic  source  polarized  with  a  300V  battery.   The 
active  face  of  the  source  was  3  in.  square.   The  signal 
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F(x,t)  =  Acos(nt-yx) 


6    =  Ay 
max     ' 


a  =  tan-1(D/L) 


- 


-! 


jt£  =  -Aysin(f2t-Yx)  =  6 


L  = 


. 


D  = 


L,  +  L0   | 
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aL 


a  =  26 


max 


D  =  2AyL 


• 


■ 
A  =  D/(2yL) 
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FIGURE  8.   THE  SLOPE-DETECTION  OPTICAL  SYSTEM 
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voltage  was  provided  by  a  signal  generator  which  could  be 
operated  in  either  CW  or  pulsed  modes.   The  CW  mode  was 
used  when  taking  data  and  the  pulsed  mode  was  used  when 
"tuning"  the  system  for  minimum  acoustic  reflection  from 
the  far  end  of  the  guide.   The  system  was  "tuned"  by  ad- 
justing the  number  and  position  of  triangular  sound-absorbent 
rubber  wedges  within  the  throat  of  the  exponential  horn 
until  minimum  reflection  was  observed. 

The  acoustic  probe  (CELESCO  Model  LC5-2)  was  a  ceramic 
hydrophone  with  a  useful  frequency  range  of  1  -  600  kHz, 
and  a  nominal  sensitivity  of  -126  dB  RE  lV/microbar  within 
this  range.   The  signal  from  the  hydrophone  was  amplified 
by  40  dB  to  provide  proper  oscilloscope  presentation.   A 
bandpass  filter  was  used  to  reject  signals  below  9  kHz  and 
above  40  kHz.   Figures  9,  10  and  11  show  the  system. 

B.   DESCRIPTION  OF  THE  EXPERIMENTS 

The  object  was  to  study  the  fluctuations  in  the  acoustic 
pressure  amplitude  of  a  traveling  wave  resulting  from  the 
interaction  of  the  acoustic  field  with  a  sinusoidal  surface 
wave  in  the  case  where  both  surface  and  acoustic  waves  were 
traveling  v/aves.   In  particular,  the  existence  of  predicted 
resonances  was  to  be  verified,  and  a  single  resonance  was  to 
be  examined  quantitatively  for  comparison  with  theory. 

The  first  experiment  was  the  investigation  of  the  resonance 
conditions.   The  relationship  between  the  velocity  V  of  the 
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FIGURE  9 
THE  SYSTEM 
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FIGURE  11.   WAVEGUIDE,  BEACH  AND  OPTICAL  SCREEN 
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surface  wave  and  the  surface  wave  number  y  [9]  is  given  by 


V2  =  (g/y  +  ^1)    tanh  yh  (3.1) 


Since  V  =  P./y, 


1  1 

2  2 

H  =  y[g/y  +  Ty/p]  [tanh  yh]   .  (3.2) 


Using  Eq.  3.2,  y  vs  f  =  9./2i\   was  computed  and  is  plotted 

5 

in  Figure  12.   Equations  2. 32 and  2.40  were  then  utilized 
in  that  order  to  compute  the  required  resonant  acoustic 
frequency  for  the  combination  of  normal  mode  and  perturba- 
tion resonance  being  investigated.   A  plot  of  k  vs  f  for 
the  various  modes  is  shown  in  Figure  13.   These  computations 
were  expedited  with  the  help  of  a  Hewlett  Packard  Model 
HP9100A  programmable  desk  calculator. 

Starting  with  the  (1,1)  normal  mode  and  a  given  y,  the 
values  of  the  surface  frequency  and  acoustic  frequency  were 
computed  for  the  lowest  perturbation  resonance,  m/n  =  1/1. 
The  surface  wave  generator  was  set  to  the  selected  frequency 
by  counting  and  averaging  the  number  of  cycles  over  a  two 
minute  period.   The  V/AVETEK  was  then  tuned  through  the  general 
vicinity  of  the  predicted  acoustical  frequency  until  a 

resonant  condition  was  observed.   The  empirical  values  of 

2 
(k  /y)  and  (k _/y)   were  then  computed  and  plotted. 
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Using  the  predicted  values  of  acoustic  frequencies  for 
the  higher  perturbation  resonances,  e.g.,  m/n  =  2/1  and 
m/n  =  3/1,  for  the  (1,1)  normal  mode,  the  process  was  repeated 

for  the  same  surface  frequency.   Five  surface  wave  frequen- 

2 
cies  corresponding  to  (k_/y)   =  1/8,  1/4,  1/2,  1,  and  2  were 

used.   These  values  were  selected  to  give  a  representative 

2 

spread  on  the  (k  /y)  vs  (k  /y)   curves. 

X  z 

The  values  of  the  acoustical  frequencies  for  the  (1,1) 
normal  mode,  m/n  =  2/1  perturbation  resonance  were  above  the 
cut-off  frequency  for  the  (1,2)  and  (2.1)  normal  modes. 
Those  for  the  (1,1)  normal  mode,  m/n  =  3/1  perturbation 
resonance,  were  above  the  cut-off  frequency  for  the  (2,2) 
normal  mode.   When  operating  at  an  acoustic  frequency  which 
was  high  enough  to  permit  propagation  of  higher  normal  modes, 
interference  from  these  higher  modes  was  observed.   Examina- 
tion of  Figs.  1*1  -  19  and  Eq.  2.16  reveals  co-Incident  plots 
of  resonant  acoustic  frequency  f  vs  surface  wavenumber 
for  several  combinations  of  normal  mode  and  perturbation 
resonance.   For  example,  normal  mode  (1,1),  perturbation 
resonance  2/1  is  degenerate  with  normal  mode  (2,1),  perturba- 
tion resonance  1/2.   The  remainder  of  these  combinations  are 
listed  in  Table  1.   It  was  necessary  to  device  means  to 
suppress  interference  between  the  degenerate  modes.   The 
acoustic  probe  was  positioned  to  take  advantage  of  the  symme- 
try in  modal  patterns  as  shown  in  Figure  20.   When  the  (1,1) 
normal  mode  was  being  investigated,  the  probe  was  placed  on 
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TABLE  1 


DEGENERATE  MODE  COMBINATIONS 


NORMAL     PERTURBATION  NORMAL      PERTURBATION 

MODE        RESONANCE  MODE         RESONANCE 


(1,1) 

2/1 

(1,2) 

1/1 

(1,2) 

2/1 

(1,1) 

3/1 

(2,1) 

3/2 

(1,3) 

1/k 

degenerate  with  (2,1)  1/2 

degenerate  with  (2,1)  2/2 

degenerate  with  (2,2)  1/2 

degenerate  with  (3,1)  1/3 

degenerate  with  (3,1)  2/3 

degenerate  with  (3,1)  1/1 


HH 


the  center  axis  of  the  water  column  which  was  a  node  for  the 
(1,2),  (2,1)  and  (2,2)  normal  modes.   The  only  degeneracy 
which  could  not  be  suppressed  in  this  manner  was  the  (1,1) 
normal  mode,  m/n  =  3/1  perturbation  resonance  which  was 
coincident  with  the  (3,1)  normal  mode,  m/n  =  1/3  perturba- 
tion resonance  as  shown  in  Figure  21.   However,  it  was  felt 
that  the  symmetry  of  the  source  combined  with  the  relative 
strength  of  these  two  modes  made  the  contribution  of  the 
(3,1)  normal  mode  insignificant. 
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FIGURE    21.       UNSUPPRESSED   MODAL    PATTERNS 
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To  permit  investigation  of  the  higher  modes,  the  source 
was  modified  to  suppress  the  (1,1)  normal  mode  and  enhance 
the  mode  in  question.   A  sheet  of  1/4  in.  neoprene  foam  was 
placed  over  the  left  half  of  the  source  to  enhance  the  (1,2) 
normal  mode,  then  over  the  upper  half  to  enhance  the  (2,1) 
normal  mode.   A  checkerboard  pattern  was  used  to  enhance 

the  (2,2)  normal  mode.   Figure  22  shows  the  data,  plotted 

2 
on  the  theoretical  (k  /y)  vs  (k  /y)   curves. 

A  second  portion  of  the  first  experiment  consisted  of 
searching  for  unexplained  resonances.   With  the  surface 
frequency  fixed,  the  input  signal  was  tuned  from  the  (1,1) 
cut-off  frequency  to  the  (2,2)  cut-off  frequency  to  see  if 
any  unexplained  resonances  could  be  located.   None  were 
found. 

The  objective  of  the  second  experiment  was  to  compare 
the  behavior  of  the  system  in  the  vicinity  of  a  resonance 
with  theoretical  predictions.   The  study  was  limited  to 
the  (1,1)  normal  mode  at  the  perturbation  resonance, 
m/n  =  1/1  to  preclude  interference  from  higher  normal  modes. 
The  surface  wave  attenuation  was  determined  by  measuring  the 
wave  height  from  HO    cm  to  120  cm  from  the  driver  wedge  In 
5  cm  increments.   The  slope,  6,  of  the  semi-log  plot  of  the 
surface  wave  amplitude  vs  distance  from  the  source  was  then 
computed.   This  process  was  repeated  at  the  five  surface 
wave  frequencies  examined  in  the  first  experiment. 

At  the  five  corresponding  acoustic  resonance  frequencies, 
the  acoustic  pressure  probe  was  traversed  over  the  same  range 
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as  the  wave  height  measurement.   The  surface  was  unperturbed, 
Since  there  was  some  reflection  of  acoustical  energy  from 
the  termination,  readings  were  taken  at  nodes  and  antinodes 
of  the  standing  wave.   The  amplitude  of  the  readings  at 
the  nodes  were  plotted  on  a  semi-log  plot  vs  the  distance 
from  the  transducer.   The  slopes,  n ,  of  these  plots  were 
then  computed  and  compared. 

One  resonance  was  selected  for  qualitative  and  quantita- 
tive analysis.   The  surface  wave  attenuation  constant  6  and 
acoustic  absorption  coefficient  in  the  x-direction  r\   were 
measured  both  before  and  after  data  were  taken.   The  absorp- 
tion coefficient  in  the  z-direction  3  was  then  computed 
using  Eq.  2.25. 

With  the  wave  generator  set  at  a  frequency  of  4.2  Hz 
the  acoustic  signal  was  set  at  various  frequencies  in  the 
vicinity  of  the  resonant  frequency  of  13-25  kHz.   The  fre- 
quency band  of  interest  was  12.5  -  1^ . 0  kHz.   The  MEMOSCOPE 
was  triggered  in  single  sweep  mode  and  amplitude  data  were 
read  from  the  CRT  display.   Figure  23  is  a  graphical  illus- 
tration of  the  display  and  Figures  24  and  25  show  typical 
examples  of  strongly  and  weakly  modulated  signals.   The 
data  taken  were  the  values  of  2a  +  2b  and  2a  -  2b.   The 
modulation  parameter,  M  =  b/a  was  then  computed  as  follows 

X  =  [1  -  H   ;  Hi   x  100  (3-3) 
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FIGURE  23.   MEASUREMENT  OF  MODULATION  PARAMETER 
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FIGURE  2k 

STRONG 
MODULATION 
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FIGURE  2  5 

WEAK 

MODULATION 
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Eight  sets  of  data  were  taken.   The  acoustic  probe  was 
positioned  at  each  of  four  distances  from  the  source:  k2 , 
67,  90,  and  100  cm.   Four  runs  were  completed  at  the  42  cm 
position,  one  at  the  67  cm  position,  two  at  the  90  cm  position, 
and  one  at  the  100  cm  position.   The  data  are  shown  in 
Figures  26  -  33. 

C.   ERROR  ANALYSIS 

The  perturbation  constant,  e,  which  is  the  ratio  of 
waveheight  to  water  column  depth  depends  upon  how  precisely 
the  waveheight  could  be  determined.   The  deflection  of  the 
laser  beam  on  the  wall  screen  could  be  read  accurately  to 
±1  cm  and  the  path  length  determined  to  ±3  cm.   Assuming 
that  the  waveguide  was  uniform  to  within  ±1  %   of  the 
measured  dimensions,  c   becomes 


.  D  ±  1 /O  EN 

eo  "  U„  ±  0.01  I    )2Y(L  ±  3) 


where  all  dimensions  are  in  centimeters. 

This  yields  an  uncertainty  in  the  measurement  of  eQ 
of  ±5  %.      The  experimental  determination  of  the  surface  wave 
attenuation  constant,  6,  resulted  in  a  ±5  %   error  from 
"worst  case"  graphical  interpretations.   The  in  situ  value 
of  e  was  dependent  upon  the  relationship 
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-<5x 
e  =  eoe   X  (3.6) 


This  resulted  in  a  cumulative  error  in  the  in  situ  measurement 
of  ±7  %. 

The  assumed  1  %    deviation  in  the  dimensions  of  the  wave- 
guide leads  to  a  ±1.5  %   error  in  the  value  of  k  .   Examination 
of  Eq.  2.13 


1  -  2k  /y 
Kz  "    kz  J1 V  (3'7) 

(k  /y) 

shows  that  an  error  of  approximately  ±3  %   in  K  could  result. 

Li 

At  resonance,  the  spatial  decay  in  the  z-direction  can 
be  written 


(2n  +  6)k 

6  =  * -  -     (3.8) 

Kz 


The  value  of  the  acoustic  absorption  coefficient  in  the 
x-direction,  r\ ,    was  of  the  order  2  X  10   .   During  the 
measurement  of  n  it  was  not  possible  to  discern  nodes  and 
antinodes  since  the  standing  wave  ratio  was  very  nearly  one. 
Since  the  slope  of  the  semi-log  plot  of  acoustic  amplitude 
vs  distance  was  so  slight  it  was  possible  to  make  a  "worst 
case"  error  estimate  of  the  value  of  n  of  ±100  %.      However, 


the  order  of  magnitude  of  the  surface  wave  attenuation  con- 

_2 
stant,  6,  was  3  X  10   ,  large  when  compared  to  n .   Thus,  the 

combined  error  of  the  term  (2n  +  <5 )  was  approximately  8  %. 
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The  equation  for  the  modulation  parameter,  M,  can  be 
written 

,  sin(7Tz/£  ) 

M  =  [~-eTTexp(-yx)] (3.9) 

d  1 

[(akz£z)2  +  (e*z)2] 

The  term  sinCirz/Jl  ),  when  expanded,  leads  to  an  approximate 
error  of  1  %    from  dimensional  considerations.   Thus,  at 
resonance,  the  total  error  in  M  was  approximately  11  %, 
considering  all  sources  of  error  to  be  random.   In  Figures 
26  -  33,  this  11  %   error  is  shown  on  the  dashed  plot  of  the 
theoretical  modulation  parameter  as  vertical  arrows.   The 
error  flags  plotted  on  representative  data  are  indicative  of 
a  12.5  %   measurement  error  at  low  modulation,  and  a  6.5  % 
error  at  high  modulation. 

There  was  a  variation  in  the  modulation  parameter  at 
the  resonant  frequency  which  could  be  attributed  to  a  standing 
surface  wave.   The  theoretical  predictions  for  a  standing 
surface  wave  were  not  developed  in  this  work.   Unfortunately 
it  had  been  decided  to  take  all  data  at  positions  of  the 
relative  maximum  modulation.   These  positions  were  seen  to 
be  separated  by  integer  multiples  of  the  half  wavelength  of 
the  surface  wave,  consistent  with  the  presence  of  a  standing 
surface  wave  in  addition  to  the  traveling  wave.   The  varia- 
tion of  the  observed  modulation  parameter  appeared  to  range 
from  an  estimated  10  %   at  the  42  cm  position  to  as  high  as 
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50  %   at  the  90  and  100  cm  positions.   Representative  data, 
corrected  for  this  variable  error  are  plotted  as  X's  in 
Figures  26  -  32. 

The  variation  in  the  maximum  modulation  parameter  was 
the  greatest  source  of  error.  While  the  error  flags  on  the 
theoretical  curves  and  data  are  indicative  of  the  exactness 
of  the  waveguide  and  preciseness  of  the  measurements  taken, 
there  was  little  basis  upon  which  an  estimation  of  how  well 
the  traveling  wave  condition  had  been  achieved. 
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IV.   RESULTS  AND  CONCLUSIONS 

The  experimental  results  consist  primarily  of  two  sets 
of  data.   The  first  set,  shown  in  Figure  22,  shows  excellent 
correlation  with  theory  and  proves  that  there  exists  an 
infinite  number  of  perturbation  resonance  conditions  for  the 
modulation  of  the  acoustic  signal.   The  fundamental  pertur- 
bation resonance  occurs  when  the  surface  wavelength  is  one- 
half  the  acoustic  wavelength.   The  higher  perturbation 
resonances  occur  for  a  linear  relationship  between  the 
dimensionless  parameters  (k  /y)  and  (k  /y)  .   Although  a 
complete  set  of  integer  values  of  m/n  was  not  investigated 
due  to  frequency  considerations,  it  can  be  implied  from 
the  data  presented  that  perturbation  resonances  occur  at 
all  integer  combinations  of  m/n. 

The  absorption  coefficient  of  the  surface  wave  increases 
linearly  with  frequency  as  shown  in  Figure  3^.   The  data, 
taken  over  a  period  of  several  days,  had  a  96  %   correlation 
to  the  straight  line  plotted  with  the  data.   At  frequencies 
lower  than  those  shown  on  the  graph,  a  sizable  standing  sur- 
face wave  developed,  making  it  difficult  to  observe  the 
absorption  coefficient.   Due  to  the  relatively  short  path 
length  for  the  acoustic  wave,  the  relationship  between  the 
acoustic  signal  in  the  x-direction  and  the  acoustic  frequency 
was  not  determined. 

There  exist  many  degenerate  perturbation  modes  where 
the  acoustic  resonant  frequency  and  the  surface  wave  number, 
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SURFACE  WAVE  FREQUENCY,   f  ,   (Hz) 

FIGURE    3^ 


6  it 


Y,  are  coincident  for  different  normal  modes.   These  degen- 
eracies are  listed  in  Table  1.   In  addition  to  these  modes, 
there  are  several  singular  degeneracies  where  curves  cross 
at  only  one  point. 

There  is  excellent  correlation  between  the  data  presented 
In  Figures  26  -  29  and  the  predicted  behavior  of  the  modula- 
tion parameter  in  the  vicinity  of  resonance.   There  is  con- 
sistency in  the  bandwidth  of  all  the  other  data  sets  and 
the  theoretical  bandwidths;  but,  as  can  be  seen  in  Figures 
30  -32  there  is  considerable  deviation  in  the  amplitude  of 
the  predicted  and  measured  modulation  at  the  90  cm  and  100  cm 
positions.   Two  factors  may  have  contributed  to  this  variation 

1.  Beyond  the  90  cm  position,  it  was  noted  that 
the  surface  wave  attenuation  constant  was  no  longer  linear, 
indicating  an  apparent  decrease  in  the  rate  of  decay  of  the 
surface  wave. 

2.  The  standing  wave  noted  previously  at  the  low 
surface  wave  frequencies  was  probably  present,  though  unde- 
tected at  the  higher  frequency  used  in  this  section  of  the 
experiment . 

It  would  be  expected  that  standing  wave  effects  would 
not  more  than  double  the  modulation  parameter  even  at  their 
worst,  whereas  a  factor  of  approximately  four  was  experienced. 
This  is  an  area  which  should  be  examined  more  fully  in  future 
experiments. 
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The  deviation  of  data  below  the  predicted  modulation 
parameter  seen  in  Figure  33  was  due  to  the  presence  of 
bubbles  on  the  acoustic  probe.   The  pressure-release  boun- 
dary created  on  the  surface  of  the  probe  by  the  bubbles  had 
the  expected  effect  of  decreasing  the  sensitivity  of  the 
probe  and  a  decrease  in  signal  amplitude  resulted. 

There  are  several  considerations  which  should  be  made 
before  next  attempting  this  experiment : 

1.  Extreme  care  should  be  taken  to  ensure  that  the 
waveguide  dimensions  are  exact.  This  could  be  accomplished 
by  constructing  the  support  walls  of  metal,  which  should 

be  machined  if  necessary  to  produce  absolutely  uniform  dimen- 
sions.  In  addition,  the  pressure  release  walls  should  be 
made  from  styrofoam,  planed  to  desired  dimensions,  instead 
of  polyurethane  which  is  not  as  sturdy. 

2.  More  careful  attention  should  be  given  to  the 
construction  of  a  termination  which  would  more  adequately 
absorb  surface  wave  energy.   The  exponential  horn  could  be 
lengthened  and  the  surface  wave  could  be  channeled  to  many 
small  absorption  chambers  by  means  of  thin,  sound-absorbent 
dividers.   It  is  believed  that  this  modification  would  not 
adversely  effect  the  excellent  acoustic  termination  achieved 
in  this  experiment. 

3.  A  design  for  automated  data  taking  should  be  con- 
ceived.  This  would  require  a  detector  to  separate  the  modula- 
tion from  the  acoustic  signal,  and  a  means  for  automatically 
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recording  the  modulation  as  the  probe  is  traversed  down  the 
waveguide. 
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